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£T) ' Abstract 

o ' 

, We present new sets of fragmentation functions for charged pions and kaons, both 

at leading and next-to-leading order. They are fitted to data on inclusive charged- 
r S^ • hadron production in e + e~ annihilation taken by TPC at PEP (i/i = 29 GeV) and 

to similar data by ALEPH at LEP, who discriminated between events with charm, 
C^h! bottom, and light-flavour fragmentation in their charged-hadron sample. In con- 

trast to our previous analysis, where we only distinguished between valence-quark, 
sea-quark, and gluon fragmentation, we are now able to treat all partons indepen- 
dently and to properly incorporate the charm and bottom thresholds. Due to the 
sizeable energy gap between PEP and LEP, we are sensitive to the scaling violation 
in the fragmentation process, which allows us to extract a value for the asymptotic 
scale parameter of QCD, A. Recent data on inclusive charged-hadron production 
in tagged three-jet events by OPAL and similar data for longitudinal electron po- 
larization by ALEPH allow us to pin down the gluon fragmentation functions. Our 
new fragmentation functions lead to an excellent description of a multitude of other 
e + e~ data on inclusive charged-hadron production, ranging from y/s = 5.2 GeV to 
LEP energy. In addition, they agree nicely with the transverse-momentum spectra 
of single charged hadrons measured by HI and ZEUS in photoproduction at the ep 
collider HERA, which represents a nontrivial check of the factorization theorem of 
the QCD-improved parton model. In this comparison, we also find first evidence 
for the interplay between the direct- and resolved-photon mechanisms and for the 
existence of a gluon density inside the photon. 



* Supported by Bundesministerium fiir Forschung und Technologie, Bonn, Germany, under Contract 
05 6 HH 93P (5), and by EEC Program Human Capital and Mobility through Network Physics at High 
Energy Colliders under Contract CHRX-CT93-0357 (DG12 COMA). 
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1 Introduction 



The advent of precise data on inclusive hadron production at pp, ep, and e + e~ colliders 
offers the challenging opportunity to test quantitatively the QCD-improved parton model. 
The inclusive cross section is expressed as a convolution of the parton density functions 
(PDF), the partonic cross sections, and the fragmentation functions (FF) of the quarks 
and gluons into the outgoing hadrons. The factorization theorem ensures that the PDF 
and FF are universal and that only the hard-scattering partonic cross sections change 
when different processes are considered. While there exist highly sophisticated sets of 
the PDF for the incoming protons and photons, which are needed for the analysis of 
inclusive particle production in pp and ep reactions, the status of the quark and gluon FF 
is much less advanced. The most direct way to obtain information on the FF is to analyse 
the energy spectrum of the hadrons into which the jets produced by e + e~ annihilation 
fragment. 

Recently, we performed such an analysis to extract FF for charged pions and kaons at 
leading order (LO) and next-to- leading order (NLO) in QCD 0. These FF were generated 
through fits to e + e~-annihilation data taken at centre-of-mass (CM) energy yfs = 29 
GeV by the TPC Collaboration at SLAC 0. About one year ago, these were the most 
precise data for charged-pion and -kaon production. Our FF also led to rather satisfactory 
descriptions of other e + e~ data on charged-particle production at lower (DORIS [§,f|), 
similar (PEP f§), and higher energies (PETRA g|7fl, TRISTAN @, and LEP @). 

These new parameterizations — in particular, the NLO sets — were tested against data 
on single-charged-hadron production obtained by the HI Collaboration in ep-scattering 
experiments at HERA [10] in Refs. 0Jll| and against similar data obtained in various pp 
experiments in Ref. The agreement between experimental data and theory turned 
out to be very satisfactory. This motivates us to further improve the FF, especially since 
detailed high-statistics data on charged-hadron production at LEP have become available 
T3| . p~^ , p~5| , [T6|1 . For our extended analysis, we shall employ very accurate data from the 



ALEPH Collaboration at LEP, who partly discriminated between pions and kaons [O] . In 
the data sample without discrimination of the hadron species, they distinguished between 
three cases, namely fragmentation of (i) u, d, and s quarks, (ii) b quarks only, and (iii) all 
five quark flavours (w, d, s, c, and b). In particular, the latter data on the fragmentation 
of specific quark flavours enables us to remove a theoretical assumption which we had to 
make in our earlier work [|IJ, namely that the s, c, and b (d, c, and b) quarks fragment 
into charged pions (kaons) in the same way. As is well known, b quarks are produced 
considerably more copiously at the Z resonance than at lower CM energies. Therefore, 
ALEPH was able to separate the fragmentation of b quarks into charged hadrons from the 
fragmentation of all other quarks. On the basis of these new ALEPH data on charged- 
hadron production, Cowan |13| has performed a similar analysis using the Altarelli-Parisi 
(AP) equations [[nj for the Q 2 evolution and incorporating data from other experiments 
at lower energy. The scope of his analysis is quite different from ours. He concentrates 
on charged-hadron data and is chiefly interested in establishing the scaling violation so as 
to determine a s (M|). Another significant difference is that he uses a rather high starting 
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scale, Qo = 22 GeV, which lies well above the scale characteristic for our ep studies. 
Similarly to the case of the PDF, the method of choice for our purposes it to use simple 
parameterizations at the low threshold scales. 

A problem that requires special attention is related to the gluon fragmentation into 
charged hadrons. Although the gluon does not couple directly to the electroweak currents, 
it contributes in higher orders and mixes with the quarks through the Q 2 evolution. This 
renders it very difficult to pin down the gluon FF by using just the energy spectrum of 
hadrons produced by e + e~ annihilation. However, the gluon fragmentation can be probed 
in e + e~ annihilation by exploiting tagged three-jet events or longitudinal polarization. As 
in our earlier work [II], we use the information on tagged three-jet events made available 



by the OPAL Collaboration |L5] to test our gluon FF. Recently, in Ref. ||12|| , it was 
shown that a sufficiently sizeable gluon FF into charged hadrons is essential to reach 
good agreement with the data of single-charged hadron production in pp processes. Since 
gluon production is also significant in ep collisions at HERA, we should be able to further 
test the strength of the gluon FF extracted from e + e~ data by making comparisons with 
data on inclusive photoproduction of charged hadrons recently collected by the HI [|T(| 
and ZEUS Collaborations 0. 

It is the purpose of this work to make use of the recent data by ALEPH [13,14| in 
addition to the data by TPC || and OPAL [15 to construct new sets of LO and NLO FF 
without imposing identities between the FF of different flavours. (We shall still identify 
the FF of u and d quarks into charged pions and those of u and s quarks into charged 
kaons.) The new FF sets will be tested against additional e + e~ data from LEP and 
e + e _ colliders with lower CM energies. In addition, we shall study single-charged hadron 
production in ep collisions under HERA conditions in order to check the consistency of 
our FF with the new HI [pij and ZEUS data |TH|, in particular with respect to the issue 
of the gluon FF. 

This work is organized as follows. In Section 2, we shall introduce the formalism 
needed to extract the FF from e + e~ data on inclusive single-hadron production at NLO. 
In Section 3, we shall present the formulae that describe inclusive single-hadron production 
in ep collisions with almost real photons. Section 4 will deal with the actual analysis and 
the discussion of our results. In Section 5, we shall check our results against e + e~ data 
which we did not use in our fits and against the recent HI and ZEUS data. Our conclusions 
will be summarized in Section 6. As in our earlier work we shall list, in the Appendix, 
simple parameterizations of our FF sets. 



2 Formalism for e + e reactions 

The inclusive production of a hadron h by e + e~ annihilation, 

e+e" ^(j,Z)^h + X, (1) 

is completely characterized by three observables. These are the energy fraction of the 
outgoing hadron x = lE^j y/s, its angle with the beam axis 9 and the total energy of the 
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system ^/s. The cross section of process ([!]) exhibits the angular structure 

d 2 a(e + e~ -> h + X) 3, 2 ^da T 3 2/1 dV L 3 do A 

v i = - 1 + CO s 2 fl — - + -sin 2 #— - + -cos#^, 2 

dx dcosO 8 dx 4 dx 4 ax 

where the superscripts T and L denote the contributions due to transverse and longitu- 
dinal polarizations, respectively, and the asymmetric term, labelled A, accounts for the 
interference of the photon with the Z boson. Usually, experiments do not determine 9 
distributions. Thus we shall integrate over 9. This eliminates the asymmetric term in 
Eq. (§. 

The partonic subprocesses may be exactly treated in perturbative QCD. However, 
this is not yet possible for the process of hadronization. As explained in the Introduction, 
the latter is described in terms of phenomenological FF, which must be extracted from 
experiment. In the language of the QCD-improved parton model, the x distribution of 
process (|l|) emerges from the x distribution, (da a /dx)(x, fi 2 , Q 2 ), of e + e~ — > a+X through 
convolution with D^(x,Q 2 ), 

1 da(e + e~ — > h + X) ^ f dz h( 2 . 1 da a ( x 



<?tot dx V J z 7 o tot dy 

Here, the sum extents over all active partons (a = u, d, s, c, b, g), fi is the renormalization 
scale of the partonic subprocess, and Mf is the so-called fragmentation scale. At NLO, 
Mf defines the point where the divergence associated with collinear radiation off parton 
a is to be subtracted. Throughout this paper, we neglect finite-quark-mass effects in the 
matrix elements. In order to be able to compare our calculations with experimental data, 
we normalize the cross section to 



N -l 

i=l 



Otot = N c ^e 2 q .a , (4) 



where a = (47ra 2 /3s) is the total cross section of e + e" — > for massless leptons, 

N c = 3, e qi is the electroweak charge of quark q± in units of the positron charge, and 
Nf is the number of active flavors. In Eq. (|4]), we take into account the running of the 
fine-structure constant, a, and the enhancement of the contributions due to down-type 
quarks at the Z pole. Since the absolute normalization is determined experimentally only 
to the 4% level, we neglect corrections to cr to t in higher orders of a s . 

To NLO in the MS scheme, the cross sections of the relevant subprocesses are given 

by 



— ^(y, A 2 , Mj) = 4 t N c — {s(l - y) + °" {p2) 



<?tot dy 1 qi a tot { 2tt 

-^p(y,^ 2 ,M 2 ) = 2 a ^ 2) 

rr . . nit * 



P£> T) (y)^ + Kj(y) + Ki(y) 



a to t dy 3 2n 



Pi°/ ) (y)^^ + Kj(y) + K^(y) 



(5) 
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Pah (v) i s the LO term of the transverse a — > b splitting function, 



Pl b {y^ s {Q 2 ))=PT\y) + 



a s (Q 2 
2tt 



(6) 



The K functions have been presented in Ref. ||19|| . The NLO formula of a s (fi 2 ) may be 
found, e.g., in Ref. [20[|. The scales \x and Mf are usually identified with the only intrinsic 
energy scale, yfs. For this choice of scales, the terms in Eq. ([5]) involving ln(s/M 2 ) are 
suppressed, so that the NLO corrections are expressed just in terms of the K functions. 
At LO, only the transversely polarized photon contributes to the cross section. The 
longitudinal cross sections of the subprocesses are thus given by 



1 da^. , 9x 
crtat dy 



C F e\N c ^ 
2vr q% a tot 

asifi 2 )^ 4(1 -y) 



2tt 



c f 



(7) 



Having defined the partonic subprocesses, we now turn to the FF. Their x distributions 
are not yet calculable in the framework of perturbative QCD. However, once we know 
them at some scale Qg, the Q 2 evolution is determined by the AP equations [jT7| . Our 
task is thus to construct a model for the x distributions at a starting scale Qq, which, 
after evolution, fits the data at scale Q 2 . The AP equations read 



d 



d\nQ- 



:D h a (x,Q 2 



a s (Q 2 
2tt 



b Jx v \y 



That is, we have to solve a system of integro-differential equations. This may be achieved 
with the help of the Mellin-transform technique |2l| . The essential property of this trans- 
formation is that it renders convolutions to products. In fact, denoting the moments of 
Pj b and D% by A ab and M a , respectively, we have 



d 



d\nQ 2 
d 

dhiQ 2 
d 



M+(n,Q 2 ) 
M s (n,Q 2 ) 



MQ 2 



M G (n,Q 2 



2tt 

a s (Q 2 ) 

2tt 
a s (Q 2 ) 



-A NS (n,a s (Q 2 ))Mr(n,Q 2 



A qq (n, a s (Q 2 ))M s (n, Q 2 ) + A gq (n, a s (Q 2 ))M G (n, Q 2 ) 



f/lu() 2- — < o T ^(n^siQ^M^Q^ + A^asiQ^MG^Q 2 )], (9) 

where A^s refers to the usual non-singlet combination of quark and antiquark splitting 
functions [pi]] , the LO expressions for A a b may be found in Refs. [pl] , p^] and the NLO 
ones in Ref. [E3|, and 



M+(n,Q 2 ) 

M E (n,Q 2 ) 
M G (n,Q 2 ) 



M qt (n,Q 2 ) + M qi (n,Q 2 



2N; 



■M s (n,g 2 ), 



i=l 



M g {n,Q 2 ). 



(10) 
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We do not include combinations M„- 



Mq i} which are asymmetric under baryon- 



number symmetry, since we do not distinguish between quarks and antiquarks. Note that 
the non-singlet terms, Mf , decouple from the gluon. The solutions of Eq. (P) factorize, 



M(n,Q 2 2 ) = E(QlQ 2 2 ,N f )M(n,Ql). 



11) 



The relevant evolution operators, E(Q\,Q"^), have been calculated in Ref. 
stance, in the non-singlet case, the NLO operator reads 



211. For in- 



O A NS 



1 + 



2n 



P\ 4 (o) ^ .(1) 
W A NS jA NS 



(12) 



where (3q = (33 — 2iVy)/3 and (3\ = (153 — 19Nf)/3. The corresponding equation for the 
coupled system exhibits a similar structure and is omitted here for ease of presentation. 
The subtraction in Eq. ( |12"D has the effect that, at Q 2 = Qq, the FF are equal to the 
ansatz also at NLO. In this respect, we differ from the approach of Ref. Q. We evolve 
the FF in three steps. We start with three quark flavours at the scale Qq and evolve up 
to the charm threshold using the evolution operators E(Qq, 4m 2 ., 3). There the charm FF 
is added to the set of FF, which, in the second step, are evolved to the bottom threshold 
employing E{4m 2 , 4m 2 ,, 4). In the final step, the FF of all five quark flavours are evolved 
to the considered CM energy Q 2 with the help of E[4m\, Q 2 ,5). 

Strictly speaking, our QCD formalism is appropriate only for the strongly produced 
pions and kaons. However, a minor fraction of the observed events are due to pions 
and kaons that are produced through weak decays of primary D and B mesons. In this 
case, it is unclear whether the usual Q 2 evolution may be applied. In our analysis, we 
assumed that this may be done. A more careful study would require that the experiments 
discriminate strongly and weakly produced pions and kaons, which was not the case for 
the data used here. 



3 Formalism for ep reactions 

In this section, we shall outline the formalism pertinent to inclusive photoproduction of 
hadrons with ep colliders, and in particular with HERA. According to present HERA 
conditions, E e = 26.7 GeV electrons collide with E p = 820 GeV protons in the laboratory 
frame, so that = 296 GeV is available in the CM frame. Among the experimentalists 
it has become customary to take the rapidity, ?/i a b, of hadrons travelling in the proton 
direction to be positive. The CM rapidity, i/cm, is related to y^b by 

vcm = yiab- \ ln ^-- ( 13 ) 

In photoproduction, the electron beam acts like a source of quasi-real photons, so that 
HERA is effectively operated as a 7p collider. The appropriate events may be discrim- 
inated from deep-inelastic-scattering events by electron tagging or anti-tagging. The 
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photon flux is well approximated by the Weizsacker- Williams formula [E4 



1 + "-'''h^tri.' 1 



2 ■ e V O 2 O 2 ■ 

mm \ "B max mm , 



(14) 



where z = E 1 /E ei Q 2 min = m 2 e z 2 /(l - z), and Q 2 m3X = 0.01 GeV 2 (0.02 GeV 2 ) for 
tagged events at HI (ZEUS). The cross section of ep — > h + X emerges from the one 
of 7p — > h + X by convolution with / 7 / e (z). By kinematics, z min < z < 1, where 
-2mm = Prexp(— j/cm)/[v^ — Prexp(?/ CM )]. I n our analysis, we shall impose 0.3 < z < 0.7 
and 0.318 < z < 0.431 to be in conformity with the HI and ZEUS event-selection criteria, 
respectively. 

It is well known that 7p — > h + X proceeds via two distinct mechanisms. The photon 
can interact either directly with the partons originating from the proton (direct photo- 
production) or via its quark and gluon content (resolved photoproduction). Both con- 
tributions are formally of the same order in the perturbative expansion. Leaving aside 
the proton PDF, G^(x, M 2 ), and the FF, D^(x, M%), which represent common factors, 
the LO cross sections are of 0(aa s ) in both cases. In the case of the resolved mecha- 
nism, this may be understood by observing that the ab — > cd cross sections, which are 
of 0(a 2 ), get dressed by photon PDF, G^x, M 2 ), whose leading terms are of the form 
a ln(M 2 /A 2 ) oc a/a s , with A being the asymptotic scale parameter of QCD. Here, a, b, c, d 
denote quarks and gluons. In fact, the two mechanisms also compete with each other nu- 
merically. Resolved photoproduction dominates at small pt and positive y^, while direct 
photoproduction wins out at large pr and negative yi^. 

The LO calculation suffers from significant theoretical uncertainties connected with 
the freedom in the choice of the renormalization scale, fi, of a s (/i 2 ) and the factorization 
scales, My, M p , M h .f\ In order to obtain reliable predictions, it is indispensable to proceed 
to NLO. We shall first consider resolved photoproduction, which is more involved. Starting 
out from the well-known LO cross section of jp — > h + X, one needs to include the 
NLO corrections, K a b^ c , to the hard-scattering cross sections, to substitute the two-loop 
formula for a s , and to endow G^, Gf, and D% with NLO evolution. This leads to 

" P h a.b.r/ X h 

, (15) 



1 

X — 
7TS 



- ^pK*, v; v 2 )5(l -w) + ^^K ab ^ c (s, v, w; /i 2 , M 2 , M 2 , M 2 ) 
v dv 2tt if 



where da® b ^ c /dv are the LO hard-scattering cross sections, v = 1 + t/s, and w = —u/(s + 
t), with s = (p a + Pb) 2 , i = (p a — Pc) 2 , and u = (p b — p c ) 2 being the Mandelstam variables 
at the parton level. The parton momenta are related to the photon, proton, and hadron 
momenta by p a = x 7 p 7 , p& = x p p p , and p c = Ph/xh- The indices a, b, c run over the gluon 
and Nf flavours of quarks and antiquarks. In this paper, we shall take Nf = 5. The 

1 The LO cross section of direct photoproduction does not depend on M 7 . 
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functions may be found in Ref. p5| for M 1 = M p . This restriction was relaxed in 



26 



K ab . 
Ref. 

The NLO cross section of direct photoproduction emerges from Eq. (|T5D by substi- 
tuting G^^-y, = ^(1 — x */)i replacing da a b^ c /dv and K a ^ c by da 7 b->c/dv and K^b^c, 
respectively, and omitting the sum over a. The K 7 b->c functions were first derived in 
Ref. |27| setting M 7 = M p = and taking the spin-average for incoming photons and 
gluons to be 1/2. In Ref. [26|1 , the scales were disentangled and the spin-average conven- 
tion was converted to the MS scheme, i.e., to be l/(n — 2), with n being the dimensionality 
of space-time. Analytic expressions for the K^b^c functions are listed in Ref. |28| . Quan- 
titative studies of inclusive photoproduction of hadrons at HERA via the direct- and 
resolved-photon mechanisms in NLO may be found in Refs. |26|j28|j29| . 



4 Results 

In general, data of inclusive hadron production at e + e~ colliders are most suitable for the 
extraction of the FF; in the case of fixed-target, collider, and ep data, the information 
on the FF is obscured by theoretical uncertainties arising from the PDF and the choice 
of factorization scales connected with the initial state. For our analysis, we select the 
data on charged-pion and -kaon production taken at energy y/s = 29 GeV by the TCP 
Collaboration at SLAC [0] and those collected at y/s = Mz by the ALEPH Collaboration 



at LEP 1 14]. In order to constrain the FF of the different quark flavours, we employ 
the ALEPH data |TB| on charged- hadron production for the three cases: (1) sum over 
all quark flavours (u, d, s, c, b); (2) sum over u,d,s quarks; and (3) b quark. These data 
combine small statistical errors with fine binning in x and are thus more constraining 
than data collected by other e + e~ experiments in the energy range from 5.2 to 55 GeV. 
For this reason, our approach is to fit exclusively to TPC and ALEPH data and to use 
the other data for cross checks. To control the gluon FF, we exploit, as in our previous 
work the information on the tagged three-jet events by OPAL |15[] presented by Nason 



and Webber ]{HS|] , who applied small corrections to the original data. However, these data 
comprise just five points and are not very constraining. For the fitting procedure, we use 
the x bins in the interval between x min = max(0.1, 2 GeV/ yfs) and 0.8 and integrate the 
theoretical functions over the bin widths as in the experimental data. The restriction at 
small x is necessary to exclude events in the non-perturbative region, where mass effects 
are important and our formalism is bound to fail. We parameterize the x dependence of 
the FF at Ql as 

D h a {x,Ql) = Nx a {l-xf. (16) 

In contrast to our previous work [|J , we now only impose the conditions +7T = D£ ++7r 
and Df ++K = D^ ++K . For all the other FF, including those of the gluon, we keep N, 
a, and (3 as independent fit parameters. This means that, together with A, which we also 
keep as a free parameter, we have a total of 31 independent fit parameters. 

The quality of the fit is measured in terms of the average Xdof f° r an selected data 
points. Our technical procedure is as follows. We consider Eq. @ at /i = Mf = y/s = 
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29 GeV and Mz as a function of the 31 parameters that determine the x dependence 
of the FF at the respective values of Qq, which we take to be Qq = \/2 GeV for u, d, s, 
Q = m (r] c ) = 2.9788 GeV for c, and Q = m(T) = 9.46037 GeV [0 for b. Using a 
multidimensional minimization algorithm [fjOfl , we search this 31-dimensional parameter 
space for the point at which the deviation of the theoretical prediction from the data 
becomes minimal. When we fit to charged-hadron data, we must include the contribu- 
tion due to protons and antiprotons. To this end, we introduce a function, f(x), which 
parameterizes the ratio of the cross sections of p + p and tt + + ir~ production as measured 



by ALEPH |4 



f(x) = 0.195 - 1.35 (x-0.35) 2 . (17) 
Thus, the charged-particle production cross section is given by 

da h++h ' . x . ^da* + +^ da K++K ~ , , 

[! + /(*)] — T— + 



doc doc doc 

To test the perturbative stability, we perform fits in LO and NLO, although we are 
primarily interested in the NLO results. 

Our results are listed below. For the average of charged pions, we obtain 

-0-842 (1 _ xl.43 



D (n,LO) 
u 


V X, Qq) 


= D^ LO \x,Ql) 


= 1.09a; 


D (n,LO) 


V X, Qq) 


= 3.47a;- L037 (1 


- a:) 3 ' 88 , 


d Ulo) 


y x i Qq) 


= 4.52 a;" ' 898 (1 


- xf m , 


D («,LO) 


V X, Qq) 


= 3.62a:- L128 (1 




D (n,LO) 


y X, Qq) 


= 6.20x~ - 351 (l 


- x) 300 



(19) 



in LO and 



x) L47 , 



D^ NLO \x,Ql) =D^ NLO \x,Ql) = 1.16x-°- 762 (l 
D^ NLO \x, Ql) = 3.96 a;- - 825 (1 - x) 4 38 , 
D^ NLO \x, Ql) = 3.69 a;- - 873 (1 - a;) 5 - 06 , 
D^' NLO \x, Ql) = 4.00 a;- 1 ' 037 (1 - a;) 7 75 , 

D^ NLO \x, Ql) = 5.11 a;- - 287 (1 - a;) 2 - 61 (20) 
in NLO. In the case of kaons, we find 

34 (1- a;) 1 - 04 , 



D (K,LO) 




Ql) 


= Df< LO \x,Ql 


) = 0.38 a; 


d \k,lo) 




Ql) 


= 1.11 a;" - 942 (1 


-x) 2 - 81 , 


D (K,LO) 




Ql) 


= 0.55x-°- 764 (l 


-x) 2 - 69 , 


D (K,LO) 




Ql) 


= 0.72 a;- ' 792 (1 


-x) 2 - 81 , 


D (K,LO) 


\ X i 


Ql) 


= 0.43 a;- - 374 (1 


-x) 2 - 69 



(21) 
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TABLE I. CM energies, types of data, numbers of data points used, and Xd.o.f. va hies 
obtained at NLO and LO for the various data samples discussed in the text. The data 
used in the fits are marked by an asterisk. 



V s [OeVJ 


Data type 


Ref. 


No. of points 


„ 2 ATT ( \ 
Xd.o.f. 111 NL0 


Xd.o.f. 111 L0 


91.2 


a (all flavours) 

a h {u, d, s) * 

a [0) 

a n * 
a K * 

a h (longitudinal) 
u^ 1 (all flavours) 


[13 
[13 
[13 
[14 
[14 
[15 
[13 




23 
23 

no 

zo 
13 
13 



10 
zz 


A O /I 

0.84 
1.13 
1.06 
1.19 
1.10 
^ .Uo 
8.65 
U.yo 


0.85 
1.41 
0.o4 

1.30 
1.57 

R OO 

o.zz 

n so 
U.oz 


55.2 


a h 


i 


7 


1.18 


1.23 


A Q 7 


(X 




Q 
O 


n o/i 


n oo 
U.yz 


35.0 


a h 


i 




22 


0.64 


1.24 


34.0 


a 71 
a K 






6 
2 


0.63 
0.66 


0.54 
0.99 


29.0 


a h 

a w * 

a K * 


i 
i 

E 


14 
18 
18 


1.55 
1.30 
0.88 


2.13 
1.49 
0.75 


9.98 




i 

i 


19 
14 


1.56 
1.00 


1.88 
1.35 


5.20 


a K 


i 
§ 


5 
3 


0.93 
1.49 


0.45 
2.31 



in LO and 



1.01 



Di K > NLO \x, Ql) = Df' NLO \x, Ql) = 0.30 oT - 978 (1 - x) 
D^ NLO \x,Q 2 ) = 1.05aT a804 (l -a;) 2 ' 54 , 
D^ NLO \x, Ql) = 0.68 or ' 753 (1 - a;) 2 ' 59 , 
Dl K > NLO \x, Ql) = 0.58 a;- - 843 (1 - x) 2 88 , 

Df> NLO \x,Ql) =0.33 x~°- 351 (l-x) - 65 ( 22) 

in NLO. Here, it is understood that the Qq values refer to the individual starting points 
given above. The Xdof values achieved for the different data sets may be seen from 
Table I. 

For the reader's convenience, we list simple parameterizations of the x and Q 2 de- 
pendences of these sets in the Appendix. We believe that such parameterizations are 
indispensable for practical purposes, especially at NLO. However, we should caution the 
reader that these parameterizations describe the evolution of the FF only approximately. 
Deviations in excess of 10% may occur for x < 0.1, in particular for the gluon. While 
this kind of accuracy is fully satisfactory for most applications, it is insufficient for the 
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comparison with the high-statistics data collected at LEP. We wish to point out that all 
Xd.o.f. values presented in this paper have been computed with the full AP-evolved results, 
which have an estimated relative error of less than 0.4%. 

Our FF exhibit good perturbative stability. In general, there is only little difference 
between our LO and NLO sets. The only appreciable difference occurs for Dg. Since 
we included in our fit high-quality data from two very different energies, namely 29 and 
91.2 GeV, we are sensitive to the running of a s and are, therefore, able to extract a value 
of A^-, appropriate to five active quark flavours. We find = 195 MeV (107 MeV) 
in NLO (LO), which corresponds to a s (M§) = 0.1154 (0.1215), in good agreement with 
the value 0.120 ± 0.008 extracted from a global fit to the Z-boson observables measured 



at LEP We find that the error on due to faked scaling effects in connection 



with the distinct flavour composition at the Z resonance is of the order of 20 MeV. We 
are quite insensitive to other sources of uncertainties in A, such as the one related to the 
relative normalization of the data. 

In this context, we also wish to mention nonperturbative power corrections to the frag- 
mentation process. We assume throughout that such corrections are small and need not 
be considered when analyzing the presently available data. The good agreement between 
our results and the various measurements nicely supports this assumption. However, 
it cannot yet be excluded theoretically that corrections to inclusive timelike processes 
(such as e + e~ annihilation) might be of order l/Q, which would correspond to a sizeable 
uncertainty in A^. 

Since we have built in the cc and bb thresholds, we have three different starting scales 
Qq. To illustrate the relative size of the FF, we have plotted them in Fig. |I]a-d as functions 
of x for Q = 10 GeV. We observe that, in the case of the pion sets, the 6-quark distribution 
is rather soft, whereas the w-quark distribution is quite hard. 

The goodness of our fits to the TPC @ and ALEPH data [[TJ[] on charged-pion and 
-kaon production may be judged from Figs. ^| and The fit to the charged- hadron data 
by ALEPH is seen in Fig. |], where we show the contributions due to (1) u,d,s, (2) b, 
and (3) all flavours. We also compare the contribution due to gluon fragmentation into 



charged hadrons with the data obtained in the three-jet analysis by OPAL The full 
(dotted) curves correspond to our NLO (LO) analysis. In the case of pion and charged- 
hadron production, our fits describe the data quite well, even for x < x min . Specifically, 
we find Xd.o.f. values of 1.28 (1.32) at NLO (LO), if we include the gluon information. 
Leaving out the gluon data, we obtain 1.06 (1.14). As may be seen from Table I, the five 
gluon data points contribute an extraordinarily high value of xi o f • 

This concludes the description of our fit procedure and the quality of the fit to the 
data used. In the next section, we shall first study to what extend our FF also account for 
other data, both below and at the Z resonance, including the longitudinal cross section 
at the Z peak, which has recently been measured |]13||. In addition, we shall confront our 



prediction for charged-particle production in ep scattering at HERA with the data by HI 
and ZEUS 1X31. 
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5 Applications 



In order to gain confidence in the validity of our new FF sets, we shall now test them 
against data which were not used in our fits. Specifically, we shall consider data on 
charged-pion and -kaon production taken by the DASP and ARGUS Collaborations 
U at DORIS and the TASSO Collaboration [0] at PETRA. Furthermore, we shall make 
comparisons with charged-hadron data collected by MARK II |J at PEP, CELLO |J 
at PETRA, AMY § at TRISTAN, and DELPHI [| at LEP. In the latter case, we 
shall assume that the percentage of the produced protons and antiprotons relative to the 
charged pions may be approximated for all energies by the function f(x) given in Eq. (|I7D , 
so that Eq. ( |I8D is valid. For each of these data samples, we list in Table I the CM energy, 
the number of data points with x min < x < 0.8, and the resulting XcL.f. values in NLO and 
LO. We see that all these data are fitted rather well, with Xd o f values of the order of unity. 
In fact, these Xdoi vahies are comparable with those we obtained for the TPC and 
ALEPH ||T3| , |l4| data, to which we actually fitted. This serves as a nontrivial consistency 
check for our procedure and tells us that a global fit to all data would not further improve 
the goodness of our FF. The NLO and LO values of Xdof are very similar, with a slight 
preference for our NLO set. The good agreement with the DASP |Q, ARGUS HQ, and 
TASSO data of charged-pion and -kaon production is also demonstrated in Figs. § and 
|3], respectively. The ARGUS data are also very well described for x < x cut . In Fig. [|, we 
present the comparisons of our theoretical results with the MARK II ||, CELLO ||, and 
AMY H measurements of charged-hadron production. We see from Figs. |^-^ that the LO 
and NLO calculations agree very well at ^ = 91.2 GeV, whereas they tend to deviate at 
smaller energies, the discrepancy being up to 30% at y/s = 5.20 GeV and medium x. 

Another method to extract information on the gluon FF is to analyze the longitudi- 
nal cross section da L jdx defined in Eq. (0). Similarly the total cross section da jdx = 
do T j dx + da L / dx defined in Eq. ([3]), which we have considered so far, da L /dx may be rep- 
resented as a convolution of the FF with the longitudinal cross sections of the individual 
subpro cesses given in Eq. (|7|), 



1 da L a s f 1 dz 



cr tot dx 2tt Jx z 



2N f , 

Otot i = \ \ x 



+ 0(a 2 s ). (23) 



Since Eq. (^) is only given to LO in a s , the a s value here is not expected to be the 
same as the one in the NLO relations Thus, from measurements of da L jdx and 

da T /dx, the gluon FF can only be extracted to LO. ALEPH |]§ and OPAL |32f have 



recently presented preliminary data on da L /dx. In Fig. [5], we evaluate Eq. ( p3|) with our 
NLO FF and two-loop a s and compare the result with the ALEPH data [|13|]. The result 
obtained (full curve) falls short of the data by a factor of two. At this point, we have to 
keep in mind that Eq. (^) is a LO prediction, and one should be prepared to allow for a 
K factor, which is typically larger than one. This K factor may be simulated by changing 
the renormalization and factorization scales. Assuming a scale of 20 GeV, we are able 
to nicely describe the data. This comparison provides us with some useful check on the 
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gluon FF at low x. A similar scale change was considered by Webber in connection with 
Monte Carlo studies of this observable ||32|| . 

In Fig. |], we compared our analysis with the experimental results on gluon frag- 
mentation into hadrons obtained by OPAL Jl5]| through their three-jet analysis with the 
corrections of Nason and Webber []TSJ applied. A more direct comparison with the origi- 
nal OPAL data is feasible. For this purpose, we calculate the ratio of the inclusive scaled 
energy distribution for gluon-tagged and normal-mixture events, -Rg.tag/n.mix, as a function 
of x. This quantity is represented in our analysis by 



D h Q (x ,Q 2 )E- =I^ 

da/dx 



-"■g.tag/n.mix ' \ I 



where h stands for the sum over the charged hadrons and the quantity in the denominator 
is defined in Eq. @. The result of this comparison is presented for LO (dashed line) and 
NLO (solid line) in Fig. ||. We see that our model gives a good account of -Rg.tag/n.mix for 
medium and large x. 

The factorization theorem guarantees that the FF characterize the hadronization phe- 
nomena in a process-independent way, i.e., the FF that have been extracted from e + e~ 
data may also be used to make predictions for other types of experiments, e.g., at 77, 
ep, and hadron colliders, fixed-target facilities, etc. In the following, we shall make 
NLO predictions for inclusive photoproduction of charged hadrons at HERA and con- 



front them with recent high-statistics data taken by HI |T(| and ZEUS |fL8|| . We shall 
focus attention on the pt spectrum of the produced hadrons, averaged over a certain 
?/i a b interval. We shall work at NLO in the MS scheme with Nj = 5 quark flavours, set 
\i = M 7 = M p = M h = £p T , and adopt the photon and proton PDF from Refs. |H| and 
34| , respectively, along with our FF for charged pions and kaons. Unless stated otherwise, 



we shall choose £ = 1. We shall evaluate a s with A^ = 158 MeV f53| . As we have seen 
in Sect. 3, HI and ZEUS apply different criteria to select the photoproduction events, 
which may be simulated by appropriate choices of the photon-energy cuts and the Q^ax 
value in Eq. ([14]). Another difference is that the y^ range presently covered by the HI 
detector is —1.5 < y^ < 1-5, while the ZEUS detector covers —1.2 < y lab < 1.4. 

In Figs. [7| and |8|, we confront the HI and ZEUS data, respectively, with the corre- 
sponding NLO calculations performed with our new FF sets for charged pions and kaons. 
The agreement is almost perfect for pt^4 GeV. For higher values of pr, the central values 
of the measurements tend to depart from our predictions in both directions, but, at the 
same time, the experimental errors become larger. We shall have to await more statistics 
in the high-p-p range until we can reach final conclusions concerning the agreement of 
experiment and theory. So far, HI and ZEUS have not detected separately charged pions 
and kaons. We would like to take this opportunity and encourage these collaborations 
to tackle such a separation, since this would allow us to test the factorization theorem 
of QCD in greater detail. In Fig. || we repeat the analyses of Figs. |7| and |] for the 
average of the positively- and negatively-charged kaons (dashed lines) and compare the 
outcome with the full results for HI and ZEUS shown in Figs. [7| and || respectively (solid 
lines). We see that, in the case of HI (ZEUS), the ratio of charged kaons to charged 
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pions ranges between 53% and 70% (54% and 71%) forpr between 1.5 GeV and 15 GeV. 
Having gained a solid amount of confidence in the validity of our FF, we may ask the 
question whether the HI and ZEUS data already show evidence for the coexistence of the 



direct- and resolved-photon mechanisms. Toward this end, we compare in Fig. [10] the full 
calculations of Figs. |7] and § with the respective direct-photon contributions. The direct- 
photon contribution exhibits a discontinuity at pt ~ 3 GeV. This is due to the onset of 
the charm FF in connection with the 7$ — > cc subprocess, which is not suppressed. In 
the case of HI (ZEUS), the direct-photon contribution amounts to between 8% and 48% 
(9% and 52%) of the total result for px between 1.5 GeV and 15 GeV. Also here, we 
need more statistics at high pt before we can claim to have detected the direct-photon 
mechanism in the inclusive photoproduction of hadrons at HERA. Of course, a precise 
experimental determination of the rapidity distribution for perturbatively high values of 
Pt would serve as a more direct approach to resolve this issue. The HI Collaboration has 



already undertaken a first step in that direction |TT1 . Another interesting issue is whether 



the gluon content of the resolved photon has been established by the data. In the case 
of HI (ZEUS), it makes up between 67% and 10% (66% and 7%) of the resolved contri- 
bution and between 61% and 5% (61% and 3%) of the total one for p T between 1.5 GeV 
and 15 GeV. Thus, it seems that, in the absence of the gluon density inside the photon, 
the theoretical prediction would significantly fall short of the data at low px- Again, the 
shape of the measured rapidity spectrum would serve as an additional discriminator. 

Finally, we investigate if our FF satisfy the momentum sum rules. Guided by the idea 
that a given outgoing parton, a, will fragment with 100% likelihood into some hadron, h, 
and that momentum is conserved during the fragmentation process, we expect that 

J2 [ 1 dxxD h a (x,Q 2 ) = l (25) 

holds for any value of Q 2 . In our analysis, the left- hand- side of Eq. ( |2^| ) should be smaller 
than unity, since we consider only pions, kaons, protons, and antiprotons, which does not 
exhaust the spectrum of all possible hadrons. We approximate the contributions due to 
neutral pions, neutral kaons, and protons/antiprotons by 

D?(x,Q 2 ) = l -Dl + ^-{x,Q% (26) 

Df +R \x,Q 2 ) = Dr +K ~(x,Q 2 ), (27) 
D^(x,Q 2 ) = f(x)D: + ^-(x,Q 2 ), (28) 



respectively, where f(x) was introduced in Eq. (|T7|) . Equation (p27|) is in good agreement 
with data at various energies P,|35|, while Eq. (^) follows on from SU(2) symmetry. We 
take the lower limit of integration in Eq. (|25|) to be 0.02. This is presumably too small for 
the low energies = 2, 4, 10 GeV, which might explain why the sum-rule values are in 
some cases larger than one at these energies, in particular for our LO fit. In Table II, we 
list the results obtained with our LO and NLO FF for Q = 5, 10, 91, and 200 GeV. 



The accuracy of our results is limited due to the uncertainty in the assumptions (|2q)-(]28 
and we estimate the error to be of the order of 10%. 
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TABLE II. Left-hand side of Eq. (g§ at NLO for Q = y/2, 5, 10, 91, and 200 GeV. The 
numbers in parentheses are evaluated with our LO set. We integrate over 0.02 < x < 0.98 
and sum over pions, kaons, protons, and antiprotons. 



a 


Q [GeV] 


V2 


4 


10 


91 


200 


u 


0.83 (1.10) 


0.86 (1.03) 


0.86 (0.98) 


0.81 (0.89) 


0.79 (0.86) 


d 


0.95 (1.08) 


0.96 (1.02) 


0.94 (0.97) 


0.88 (0.88) 


0.86 (0.85) 


s 


1.08 (1.70) 


1.05 (1.52) 


1.02 (1.42) 


0.94 (1.23) 


0.91 (1.18) 


c 




0.94 (1.14) 


0.92 (1.08) 


0.86 (0.96) 


0.83 (0.92) 


b 






0.92 (1.11) 


0.85 (0.98) 


0.83 (0.94) 


9 


0.96 (0.93) 


0.96 (0.87) 


0.91 (0.82) 


0.79 (0.70) 


0.75 (0.67) 



6 Summary and Conclusions 

We presented new sets of FF for charged pions and kaons, both at LO and NLO. They 
were fitted to data on inclusive charged-pion and -kaon production in e + e~ annihilation 
taken by TPC § at PEP (^ = 29 GeV) and by ALEPH 0ITJ at LEP. In their 



charged-hadron sample ALEPH discriminated between events with charm, bottom, 
and light-flavour fragmentation. This enabled us to treat all partons independently and to 
properly incorporate the charm and bottom thresholds. We just imposed the conditions 
+7T = D^ ++7T and D^ ++K = D^ ++K , which are quite natural, since we treated the 
u, d, and s quarks as massless, so that they equally participate in the strong interaction. 
This constitutes a considerable improvement of our previous analysis Jj] based just on the 
TPC data, where we could only distinguish between valence-quark, sea-quark, and gluon 
fragmentation and neglected the heavy-flavour thresholds. In order to control the gluon 
FF, we exploited recent data on inclusive charged-hadron production in tagged three-jet 
events by OPAL [15) and similar data for longitudinal electron polarization by ALEPH 



i~3fl . Due to the sizeable energy gap between PEP and LEP, we were also sensitive to 



the scaling violation in the fragmentation process. This allowed us to extract a value for 
Ag^. We found A^ = 195 MeV (107 MeV) in NLO (LO), in good agreement with the 
global analysis of LEP data |RT . 



Although our FF were only fitted to the TPC @ and ALEPH data [pjpf, it turned 
out that they lead to an excellent description of a wealth of other e + e~ data on inclusive 
charged-hadron production, ranging from yfs = 5.2 GeV to LEP energy |||],[^,]6||^|||| . In 
fact, we always obtained Xd.o.f. values of order unity. The only exception, with Xd.o.f. ~ 



9, occurred in the longitudinal cross section [ 13 1 , which is of 0(a s ) and suffers from a 
considerable scale ambiguity. 

In order to test quantitatively the factorization theorem of fragmentation in the QCD- 
improved parton model at the quantum level, we made NLO predictions for the pt spec- 
trum of charged hadrons produced inclusively in the scattering of quasi-real photons on 
protons under HERA conditions, and confronted them with new high-statistics data by 
HI [|K]] and ZEUS [H|. Also here, the outcome of the comparison was very encouraging, 



15 



especially in the lower px range, where the experimental errors are smallest. We found 
first evidence for the interplay between the direct- and resolved-photon mechanisms in 
the inclusive photoproduction of hadrons. Our comparisons also support the notion that, 
during a fraction of time, the photon can interact with the partons inside the proton like 
a gluon. 
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A Parameterizations 

For the reader's convenience, we shall present here simple parameterizations of the x and 
Q 2 dependence of our FF.0 As usual, we introduce the scaling variable 

ln(«/A 2 )- W 
For A we use the MS value appropriate to Nf = 5 flavours, since the parameterization 

(5) 

would not benefit from the incorporation of discontinuities in s. A^ is determined from 
our fit to be 107 MeV (195 MeV) in LO (NLO). Similarly to Eqs. we use three 

different values for Qo, namely 

{V2 GeV, if a = u, d, s, g 

m(r] c ) = 2.9788 GeV, if a = c . (30) 

m(T) = 9.46037 GeV, if a = b 

This leads to three different definitions of s. For definiteness, we use the symbol s c for 
charm and Sb for bottom along with s for the residual partons. 

We parameterize our FF by simple functions in x with coefficients which we write as 
polynomials in s, s c , and Sb- We find that the template 

D(x,Q 2 ) = Nx a {l -xf (31) 

is sufficiently flexible, except for D^' 10 ^, d^' nlo ^ anc l jj(k,nlo)^ w h ere we include an 
additional factor (1 + j/x) on the right-hand side of Eq. (|3lf) . For s = s c = Sb = 0, the 
parameterizations agree with the respective ansatze in Eqs. ([TP|)-(P2]). The charm and 
bottom parameterizations must be put to zero by hand for s c < and Sb < 0, respectively. 

2 A FORTRAN subroutine that returns the FF for given x and Q 2 may be obtained from the authors 
via e-mail (binnewie@ipsl07.desy.de, kniehl@vms.mppmu.mpg.de). 
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We list below the parameters to be inserted in Eq. ([31]) for charged pions and kaons 
both at LO and NLO. The resulting parameterizations correctly describe the AP evolution 
up to 10% for Q < Q < 150 GeV and x > 0.1. Deviations in excess of 10% may occur 
for x < 0.1. 

1. LO FF for (7T+ + 7T-): 

. D^°\x,Qi) = D^ LO \x,Q*y. 

N = 1.090 - 1.085s + 0.700s 2 - 0.217s 3 
a = -0.842 - 0.920s + 0.532s 2 - 0.148s 3 
P = 1.430 + 0.434s + 0.035s 2 + 0.060s 3 (32) 

• D^ LO \x,Q 2 ): 

N = 3.470 - 2.340s + 0.309s 2 + 0.120s 3 
a = -1.037 - 0.301s - 0.038s 2 + 0.031s 3 
P = 3.880 + 0.866s - 0.137s 2 + 0.058s 3 (33) 

• D<?> LO \x,Q 2 ): 

N = 4.520 - 3.966s c + 1.453s 2 - 0.261s 3 
a = -0.898 - 0.382s c + 0.027s 3 

P = 4.690 + 0.742s c (34) 

. D^ LO \x,Q 2 ): 

N = 3.620 - 9.669s 6 + 13.160s 2 , - 5.406s^ 
a = -1.128 - 1.613s 6 + 2.715sg 
P = 7.120 - 1.081s 6 + 0.936s^ + 1.678sf 
7 = -0.164s b + 0.208s 2 (35) 

• Df> LO \x,Q 2 ): 

N = 6.200 - 13.100s + 10.940s 2 - 3.395s 3 
a = -0.351 - 1.542s + 0.319s 2 + 0.074s 3 
P = 3.000 + 1.567s - 0.725s 2 + 0.465s 3 (36) 

2. NLO FF for (vr+ +7r~): 

. D£> NLO Kx,Q 2 )=D^ NLO \x,Q 2 ): 

N = 1.160 - 1.292s + 0.975s 2 - 0.341s 3 
a = -0.762 - 1.392s + 1.151s 2 - 0.429s 3 
p = 1.470 + 0.752s - 0.350s 2 + 0.175s 3 (37) 
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D^ NLO \x,Q 2 ): 



D^ NLO \x,Q 2 ): 



D^ NLO \x,Q 2 ) 



N = 3.960 - 3.148s + 1.198s 2 - 0.267s 3 
a = -0.825 - 0.552s + 0.206s 2 - 0.070s 3 
(3 = 4.380 + 0.940s - 0.245s 2 + 0.082s 3 



N = 3.690 - 3.285s c + 1.383s 2 - 0.305s 3 
a = -0.873 - 0.481s c + 0.085s 2 
= 5.060 + 0.685s c - 0.032s 2 + 0.038s 3 



N = 4.000 - 12.965s 6 + 18.913s 2 , - 10.853s 
a = -1.037 - 1.705s 6 - 0.865s^ + 1.423sf 
= 7.750 - 2.306s 6 + 1.877sg + 0.044s 3 , 
7 = -0.116s 6 



D^ NLO \x,Q 2 ): 



N — 5.110 — 9.419s + 7.144s 2 - 2.069s 3 
a = -0.287 - 2.011s 



1.211s 2 - 0.377s 3 



P = 2.610 + 2.336s - 1.580s 2 + 0.756s 3 



LO FF for (K+ + K~ 
• Di K ' LO \x,Q 2 ) = 



Df> LO \x,Q 2 ): 



Df> L °\x,Q 2 )- 



N 
a 

(3 



N 
a 



0.380 - 0.123s - 0.015s 2 + 0.010s 3 
-1.234 - 0.106s - 0.035s 2 + 0.011s 5 



1.040 + 0.708s - 0.086s 2 + 0.048s 



>^3 



1.110 - 0.674s + 0.095s 2 + 0.028s 3 
-0.942 - 0.214s - 0.067s 2 + 0.041s 3 
2.810 + 0.852s - 0.172s 2 + 0.090s 3 



D( K > LO \x,Q 2 ): 



N = 0.550 - 0.392s c + 0.088s; 



-0.764 - 0.378s f . 



a 

(3 = 2.690 + 0.713s c - 0.103s 2 + 0.064s 3 
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0.720 - 0.495s 6 + 0.205s 3 , 

-0.792 - 0.296s 6 - 0.261st + 0.421st 

2.810 + 0.686s 6 



0.430 - 0.915s + 0.885s 2 - 0.316s 3 
-0.374 - 2.368s + 1.475s 2 - 0.317s 3 
2.690 + 1.137s - 0.934s 2 + 0.698s 3 



• Di K > NL °\x,Q 2 ) = Df' NLO \x,Q 2 ): 

N = 0.300 + 0.039s - 0.189s 2 + 0.070s 3 
a = -0.978 - 0.304s - 0.004s 2 - 0.004s 3 
P = 1.010 + 1.193s - 0.654s 2 + 0.223s 3 

N = 1.050 - 0.546s + 0.071s 2 + 0.009s 3 
a = -0.804 - 0.347s + 0.023s 2 + 0.007s 3 
P = 2.540 + 0.919s - 0.321s 2 + 0.138s 3 

• Df> NL °\x,Q 2 ): 

N = 0.680 - 0.415s c + 0.092s 2 
a = -0.753 - 0.484s c + 0.099s 2 
P = 2.590 + 0.661s c 

. D { b K,NLO \x,Q 2 ): 

N = 0.580 - 0.381s fe + 0.114s 3 
a = -0.843 - 0.414s 6 - 0.120s 2 . + 0.242s^ 
P = 2.880 + 0.614s fe 

• Df> NL °\x,Q 2 ): 

N = 0.330 - 0.196s - 0.056s 2 + 0.057s 3 
a = -0.351 - 0.026s - 0.480s 2 + 0.148s 3 
P = 0.650 + 2.326s - 0.840s 2 + 0.314s 3 
7 = 1.033s - 0.479s 2 



. Di K > LO \x,Q 2 ): 

N — 
a = 
P = 

• Df' LO \x,Q 2 ): 

N = 
a = 
P = 

NLO FF for (K + + K~): 
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FIGURE CAPTIONS 



Figure 1: x dependence of the FF at Q 2 = 100 GeV 2 for (a) charged pions at NLO, (b) 
charged pions at LO, (c) charged kaons at NLO, and (d) charged kaons at LO. 

Figure 2: Differential cross section of inclusive charged-pion production at LO (dashed 
lines) and NLO (solid lines) as a function of x at y/s = 5.2 g, 9.98 g, 29 @, 34 [|, and 



91.2 GeV [T3|. Upper curves correspond to lower energies. The theoretical calculations 
are compared with the respective experimental data. 

Figure 3: Same as in Fig. ^ for charged kaons. 

Figure 4: Differential cross section of inclusive charged-hadron production at LO (dashed 
lines) and NLO (solid lines) as a function of x at s/s = 29 @, 35 @, 55.2 g, and 91.2 GeV 
|T3| . |T5fl . Upper curves correspond to lower energies. The downmost four curves all refer 
to LEP energy. The second and third of them correspond to the contributions due to 



the u,d,s quarks and the b quark JIBf , respectively. The fourth spectrum shows the x 



dependence of the gluon FF as measured by OPAL [pL5| . The theoretical calculations are 
compared with the respective experimental data. 

Figure 5: Differential longitudinal cross section of inclusive charged-hadron production 
at NLO as a function of x at \fs = 91.2 GeV. The dashed line is obtained by choosing 
the renormalization and fragmentation scales to be 20 GeV. The theoretical calculations 
are compared with the experimental data [T3| . 



Figure 6: Ratio of the inclusive scaled energy distribution for gluon-tagged and normal- 
mixture events at LO (dashed line) and NLO (solid line) as a function of x for = 
91.2 GeV. The theoretical calculations are compared with the experimental data ||15|| . 

Figure 7: The pt spectrum of inclusive charged-hadron production as measured by ZEUS 
is compared with the NLO calculation in the MS scheme with Nf = 5 flavours using 



the photon and proton PDF of Refs. and P4] , respectively. The dashed/solid/dotted 
curves correspond to the choices £ = 0.5/1/2. 

Figure 8: Same as in Fig. ^ for HI |fL0"| . 

Figure 9: Comparison of the charged-kaon contributions (dashed lines) with the full 
results (solid lines) for the theoretical calculations in Figs. [7] and [8[ For better separation, 
the charged-kaon contributions are divided by a factor of two. The upper (lower) curves 
correspond to HI (ZEUS) conditions. 

Figure 10: Comparison of the direct-photon contributions (dashed lines) with the full 
results (solid lines) for the theoretical calculations in Figs. [7] and The steps in the 
direct-photon contributions are due to the opening of the — > cc subprocess at the 
charm-quark threshold. The upper (lower) curves correspond to HI (ZEUS) conditions. 
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